Nucleon-Nucleon Scattering From Fully-Dynamical Lattice QCD 
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One of the ultimate goals of nuclear physics is to com- 
pute the properties and interactions of nuclei directly 
from Quantum Chromodynamics (QCD), the underlying 
theory of the strong interactions. Achievement of this 
goal would reveal how nuclear processes depend upon 
the fundamental constants of nature, and would enable 
the computation of strong-interaction processes of impor- 
tance in environments not attainable in the laboratory, 
such as in the interior of neutron stars. A lone, pio- 
neering study of nucleon-nucleon (NN) scattering with 
lattice QCD was performed more than a decade ago by 
Fukugita ei o/ Q . These calculations were quenched and 
at relatively large pion masses, ttZtt^ 550 MeV. 

Lattice QCD is presently the only known method of 
calculating low-energy strong-interaction processes from 
QCD. However, it is unlikely that lattice QCD will ever 
be used to calculate the properties or interactions of nu- 
clei beyond the lightest few. In order to compute the 
properties of larger nuclei, lattice QCD calculations of 
the lightest nuclei will be performed, and then matched 
to calculations of the larger nuclei using many-body tech- 
niques such as Greens Function Monte-Carlo (GFMC), 
e.g. Ref. 0, or the No-Core Shell Model (NCSM), e.g. 
Ref. 01, including up to four-body (and possibly higher) 
interactions consistent with chiral symmetry and power- 
counting. These many-body methods have had great suc- 
cess in reproducing the properties of the light nuclei using 
high-precision phenomenological potentials as input and 
recently with the chiral potentials derived from effective 
field theory (EFT) expansions. 

In this letter we present results of the first fully- 
dynamical lattice QCD calculation of the NN scattering 
lengths in both the ^S'o-channel and — ^Di-coupled- 
channels at pion masses of rriTr ~ 350 MeV , 490 MeV and 
590 MeV in the isospin- limit. Our lattice calculations are 
performed with domain- wall 0, |^ valence quarks on the 
20'^ X 64 MILC gauge-field configurations with a lattice 
spacing of 6 ~ 0.125 fm and a spatial extent oi L ^ 
2.5 fm. The dependence of the NN scattering lengths 



upon the light-quark masses has been determined to var- 
ious non-trivial orders in the EFT expansion 0,0)01) and 
is estimated to be valid up to TOtt ~ 350 MeV. We use 
the results of our lattice QCD calculation at the light- 
est pion mass and the experimentally-determined scat- 
tering lengths at the physical value of the pion mass to 
constrain the chiral dependence of the scattering lengths 
from ttItt ~ 350 MeV down to the chiral limit. 

The NN scattering lengths are determined by comput- 
ing the energy shift of the lowest-lying two-nuclcon state 
in the finite lattice volume. To extract p cot S{p) — where 
6{p) is the phase shift — the magnitude of the center-of- 
mass momentum, p, is extracted from this energy shift 
and inserted into Liischer's relation 



pcotd{p) = 
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which is valid below the inelastic threshold. The sum in 
eq. is over all triplets of integers .i such that |j| < A 
and the limit A ^ oo is implicit [lol flTl | . The s-wave 
scattering lengths in the NN systems follow from the 
effective range expansion (ERE) 
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where a is the scattering length (using the nuclear physics 
sign convention) and r is the effective range which is typ- 
ically of order the range of the interaction, ~l/m^. 

Following the LHP collaboration (LHPC) [ialli|. our 
computation is a hybrid lattice QCD calculation us- 
ing domain-wall valence quarks from a smeared-source 
on three sets of Nf — 2 -f- 1 asqtad-improved 0, 
Hsf MILC configurations generated with staggered sea 
quarks 0|. In the generation of the MILC configura- 
tions, the strange-quark mass was fixed near its physi- 
cal value, b-nis — 0.050, (where b is the lattice spacing) 
determined by the mass of hadrons containing strange 
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quarks. The two light quarks in the three sets of configu- 
rations are degenerate (isospin-symmetric) , with masses 
bmi = 0.010, 0.020 and 0.030. Some of the domain- 
wall valence propagators were previously generated by 
LHPC on each of these sets of lattices. The domain-wall 
height is m — 1.7 and the extent of the extra dimen- 
sion is L5 — 16. The parameters used to generate the 
light-quark propagators have been "tuned" so that the 
mass of the pion computed with the domain-wall propa- 
gators is equal (to few-percent precision) to that of the 
lightest staggered pion computed with the same parame- 
ters as the gauge configurations 0|. The MILC lattices 
were HYP-blocked [l2j and Dirichlet boundary condi- 
tions were used to reduce the time extent of the MILC 
lattices from 64 to 32 time-slices in order to save time in 
propagator generation. Various parts of the lattice were 
employed to generate multiple sets of propagators on each 
lattice. We analyzed one set of correlation functions on 
564 lattices with bmi = 0.030, three sets of correlation 
functions on 486 lattices with bmi = 0.020, and one set 
of correlation functions on 658 lattices with bnii = 0.010. 
The lattice calculations were performed with the Chroma 
software suite [l^ Il9| on the high-performance comput- 
ing systems at the Jefferson Laboratory (JLab). 

The contractions necessary to produce the required 
correlation functions were performed in two stages. First, 
three propagators were contracted together at a sink with 
the quantum numbers of the nucleon, and Fourier trans- 
formed on each time slice to produce a non-lattice object 
with multiple Dirac and isospin indices (on the initial 
time slice). Second, two of these objects were contracted 
together to produce the two-nucleon correlation functions 
as a function of time-slice. For an arbitrary nucleus (or 
bound and continuum nucleons), of atomic number A 
and charge Z, there are {A + Z)l (2 A — Z)\ contractions 
that must be performed to produce the nuclear correla- 
tion function. Therefore, in the ^Sq channel there are 
48 contractions, while in the coupled channels 

system there are 36. The cleanest quantity from which 
to extract the energy-difference between the two-nucleon 
statc(s) and the mass of two noninteracting nucleons was 
found to be the ratio of correlation functions 



(3) 



where / denotes the isospin of the NN system and S 
denotes its spin. The single-nucleon correlator is 



CN{t) = 5](^(i,x) A^t(o,o)) , 



(4) 



and the two-nucleon correlator that projects onto the s- 
wave state in the continuum limit is 



J2{Nnt,^)N^{t,y)N^H0,O)Nl'H0,O)) , (5) 



where a, (3, cr, p are isospin-indices and fc, I are Dirac- 
indices. The tensor X^^pap elements that pro- 
duce the correct spin-isospin quantum numbers of two- 
nucleons in an s-wave. The summation over x (and 
y) corresponds to summing over all the spatial lattice 
sites, thereby projecting onto the momentum p = 
state. The interpolating field for the proton is pi(i,x) — 
^abcufitjx) (-u''^(i, x)C75d'^(t, x)) , and similarly for the 
neutron. The ratio of correlation functions that we ob- 
tain in the^5o channel and the '^Si —^Di coupled channels, 
at the three different pion masses, are shown in fig.^and 
fig. 12 respectively. 
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FIG. 1: The logarithm of the ratio of correlation functions 
in the channel (G^") as a function of time-slice. Each has 
been off-set vertically for display purposes. 
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FIG. 2: The logarithm of the ratio of correlation functions 
in the '^Si —^D\ coupled-channels {G"^) as a function of time- 
slice. Each has been off-set vertically for display purposes. 

The results of our calculations are shown in Table 
The scale is set via the quark- mass dependence of /tt, 
which gives b = 0.127±0.001 fm ^ (consistent with 
the Sommer scale-setting procedure used by MILC T^). 
At the pion masses used in these calculations the NN 
scattering lengths are found to be of natural size in both 
channels, and are much smaller than the L ~ 2.5 fm lat- 
tice spatial extent. It is noteworthy that our scattering 
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(MeV) aCSo) (fm) a^Si) (fm) 

353.7 ± 2.1 0.63 ± 0.50 (5-10) 0.63 ± 0.74 (5-9) 

492.5 ± 1.1 0.65 ± 0.18 (6-9) 0.41 ± 0.28 (6-9) 

593.0 ±1.6 0.0 ±0.5 (7-12) -0.2 ± 1.3 (7-12) 



TABLE I: Scattering lengths in the ^5*0 channel and in the 
^Si —^Di coupled channels. The uncertainty is statistical and 
the fitting ranges are in parentheses. There is a systematic 
error of 0.1 fm on each scattering length associated with the 
truncation of the effective range expansion; i.e. the numbers 
exhibited are for — 1/pcot 5 at the measured energy-splitting. 



physical m^n-), at a renormalization scale A ~ 350 MeV. 
In both channels the lightest lattice datum constrains the 
chiral extrapolation to two distinct bands which are sen- 
sitive to both the quark mass dependence of qa and the 
sign of the D2 coefficient. As the lattice point used to 
constrain the EFT is at the upper limits of applicabil- 
ity of the EFT, we expect non-negligible corrections to 
these regions from higher orders in the EFT expansion. 
It is clear from fig. |21 and fig. 0| that even a qualitative 
understanding of the chiral limit will require lattice cal- 
culations at lighter quark masses. 



lengths at the heaviest pion mass are not inconsistent 
with the lightest-mass quenched values of Ref. . How- 
ever, one should keep in mind the effects of quenching on 
the infrared properties of the theory [2H . 

The lowest pion mass at which we have calculated is 
at the upper limit of where we expect the EFT describ- 
ing NN interactions to be valid Sin mil mill. 

While some controversy remains regarding the details 
of the NN EFT, in our present analysis, we have con- 
strained the chiral extrapolation using BBSvK power- 
counting 113 (=KSW power-counting HI [111) and W 
power-counting jij 0| in the ^S'o-channel and BB- 
SvK power-counting in the '^Si —^Di coupled channels. 
The recent lattice QCD determinations of the light-quark 
axial-matrix element in the nucleon by LHPC ^2^] and 
its physical value are used to constrain the chiral expan- 
sion of gA- Our lattice calculations of the nucleon mass 
and pion decay constant j20| — as well as their physi- 
cal values — are used to constrain their respective chi- 
ral expansions. In addition to the quark-mass depen- 
dence these three quantities contribute to the NN sys- 
tems, there is dependence on the quark masses at next- 
to-leading order (NLO) from pion exchange, and from 
local four-nucleon operators that involve a single inser- 
tion of the light-cruark mass matrix, described by the 
"Z52" coefficients IE 0, 13 • The results of this lattice 
QCD calculation constrain the range of allowed values 
for the and consequently the scattering lengths in 

the region between ~ 350 MeV and the chiral limit, 
as shown in fig. |31 and fig. 01 With only one lattice point 
at the edge of the regime of applicability of the EFT, a 
prediction for the scattering lengths at the physical pion 
mass is not possible: the experimental values of the scat- 
tering lengths are still required for an extrapolation to 
the chiral limit and naive dimensional analysis (NDA) 
is still required to select only those operator coefficients 
that are consistent with perturbation theory. The regions 
plotted in the figures correspond to values of Co - the 
coefficient of the leading-order quark-mass independent 
local operator - and D2 that fit the lattice datum and 
the physical value, and are consistent with NDA; indeed 
we have £'2(A)m^/Co(A) ~ ±0.10 in both channels (at 
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FIG. 3: Allowed regions for the scattering length in the ^So 
channel as a function of the pion mass. The experimental 
value of the scattering length and NDA have been used to 
constrain the extrapolation in both BBSvK [25I |2^ W\ and 
W power-countings at NLO. 
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FIG. 4: Allowed regions for the scattering length in the 
'^S\ —^Di coupled-channels as a function of the pion mass. 
The experimental value of the scattering length and NDA 
have been used to constrain the extrapolation in BBSvK [2^ 
power-counting at NLO. (W counting gives essentially iden- 
tical results.) 

Without the resources to perform similar lattice QCD 
calculations in different volumes, and observing that 
most energy-splitting are positive, we have assumed that 
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we have observed scattering states in each case. Calcu- 
lations in larger volumes will be done in the future to 
verify the expected power-law dependence upon volume 
that scattering states exhibit. In addition to discriminat- 
ing between bound and continuum states, calculations in 
a larger volume would reduce the energy of the lowest- 
lying continuum lattice states, and thus reduce the un- 
certainty in the scattering length due to truncation of the 
ERE. Further improvement would result from measuring 
the energy of the first excited state on the lattice, either 
with a single source or by using the Liischer- Wolff |29| 
method. 

There are many aspects of this calculation that should 
be refined in the future. The statistics should be im- 
proved by at least an order of magnitude to have a pre- 
cise extraction of the scattering lengths from each of these 
lattices. The lattice spacing effects in the present calcu- 
lation appear at ^ 0{b^) (or exponentially suppressed 
0{b) effects), and are expected to be small. However, 
the finite lattice spacing effects should be determined by 
performing the same calculation on lattices with a finer 
lattice spacing. While it would be useful to perform this 
calculation with configurations generated from domain- 
wall quarks in the sea sector as well, the resources to 
do so are not currently available to us. However, a re- 
cent theoretical investigation [Bflj of the impact of using 
a mixed-action to compute tttt scattering has shown 
it to be small. An analogous theoretical investigation 
of mixed-action EFT for NN scattering which would 
allow a continuum extrapolation remains to be carried 
through. In addition to more precise lattice QCD calcu- 
lations through an increase in computing resources, for- 
mal developments are also required. In order to have 
a more precise chiral extrapolation, calculations in the 
various relevant EFT's must be performed beyond NLO. 
Furthermore, it is clear that lattice calculations at lower 
pion masses are essential for the extrapolation to the chi- 
ral limit, and will ultimately allow for a "prediction" of 
the physical scattering lengths. 

To summarize, we have performed the first QCD cal- 
culations of nucleon-nucleon scattering by using fully- 
dynamical, mixed action lattice QCD. This work opens 
up an unexplored area of nuclear physics as it is now 
possible to perform lattice QCD calculations of simple 
nuclear systems at pion masses within the range of valid- 
ity of the NN EFT's. 
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